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ABSTRACT 



The associations of the angular positions of background quasars with 
foreground galaxies, clusters of galaxies and quasars are often attributed to the 
statistical lensing by gravitational potentials of the matter along the lines of 
sight, although it has been known that none of the individual objects (galaxies, 
clusters or quasars) are able to fully explain the reported amplitudes of the 
quasar number enhancements. This probably arises from the fact that the 
gravitational lensing effect by the environmental matter surrounding these 
objects has been ignored. In this paper we conduct an extensive study of the 
influence of the environmental matter on the prediction of quasar enhancement 
factor by employing the spatial two-point correlation function. Assuming a 
singular isothermal sphere for mass density profile in galaxy and cluster of 
galaxies, we estimate the average surface mass density £ around galaxies, 
clusters and quasars from the galaxy-galaxy, cluster-cluster, cluster-galaxy and 
quasar-galaxy correlations. Our results show that the E induced quasar number 
enhancement in the scenario of gravitational magnification depends critically 
on the mass density parameters of galaxies (Q g ) and clusters of galaxies (Q c ) in 
the universe. For a flat cosmological model of Q = 1 the environmental matter 
can indeed play an important role in the lensing origin of the quasar-quasar 
and quasar-galaxy associations if Q g ~ Q c ~ fl , while it is unlikely that £ is 
sufficient to account for the reported quasar overdensity behind quasars/galaxies 
if galaxies and clusters of galaxies contribute no more than 25% to the matter 
of the universe. Nonetheless, the recently observed quasar-cluster associations 
on scale of ~ 10 arcminutes cannot be the result of gravitational lensing by the 
cluster environmental matter 
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1. Introduction 

Detection of the associations of the angular positions of background quasars with 
foreground galaxies and clusters of galaxies on scale of ~ 10 arcminutes in recent years 
(Fugmann 1988, 1990; Bartelmann & Schneider 1993b, 1994; Rodrigues- Williams & Hogan 
1994; Wu & Han 1995; Rodrigues- Williams & Hawkins 1995; Seitz & Schneider 1995) has 
led several studies on the statistical lensing of quasars by large-scale structures traced by 
galaxies and clusters of galaxies (e.g. Bartelmann & Schneider 1993a, b; Wu & Fang 1996). 
Indeed, gravitational matter of galaxies and clusters of galaxies alone is far from providing 
the observed amplitude of quasar overdensity to that large angular scale [see Wu (1996) for 
a recent review]. Even on small scale of a few arcseconds, the quasar-galaxy associations 
cannot be well accounted for unless an unreasonably large velocity dispersion for galaxies 
or an additional matter contribution is assumed (e.g. Narayan 1989; Schneider 1989; Wu, 
Zhu & Fang 1996). In particular, the quasar environmental matter must be invoked in the 
scenario of gravitational lensing for the explanation of the existence of four quasar pairs 
within 5" but with different redshifts among the 1000 - 2000 surveyed quasars (Burbidge, 
Hoyle & Schneider 1996). 

The spatial two-point correlation function is a quantitative description of the 
environmental matter distribution. Anderson & Alcock (1986) addressed the question 
whether clustering of galaxies described by the two-point correlation function alters the 
statistical properties of gravitational lensing. From a number of Monte-Carlo simulations, 
they concluded that the effect is too small to be significant for statistics of the multiply- 
imaged quasars in the case of galaxies acting as lenses. However, based on the N-body 
simulations of formation of galaxies, Bartelmann & Schneider (1993a) did find a correlation 
of high redshift quasars with low redshift galaxies, which stems from the magnification bias 
by galaxies and their surrounding large-scale structures. The results from these two studies 
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are not inconsistent since they investigated two different phenomena though the deflectors 
are the same (galaxies). 

It is necessary to further explore a more general question: Is gravitational lensing by 
the environmental matter able to explain the reported associations of background quasars 
with foreground objects or how large is the environmental effect on the prediction of the 
background quasar overdensity around foreground objects ? We intend to answer the 
question in this paper by an extensive analytic investigation of the magnification bias 
from the two-point correlation functions and will present elsewhere (Wu, Fang & Jing, in 
preparation) a numerical study of the issue from the N-body simulations of formation of 
clustering of galaxies and clusters of galaxies. In section 2 we estimate the mean surface 
matter density superposed statistically on galaxies, clusters of galaxies and quasars from 
the two-point correlation functions of galaxy-galaxy, cluster-cluster, cluster-galaxy and 
quasar-galaxy. We then evaluate the quasar enhancement factor near the positions of 
foreground galaxies, clusters and quasars by the environmental effect in section 3. A brief 
discussion is presented in section 4. Throughout the paper we adopt a flat cosmological 
model of Q = 1 and a Hubble constant of H = 100 h km s" 1 Mpc -1 . 

2. Environmental matter contribution from the two-point correlation function 

2.1. General consideration 

The spatial two-point correlation function can be generally written as 



at the range of separation r min < r < r max , where r is the correlation length. ^A_e( r ) 
provides a conditional probability of finding neighbor objects B (quasars, galaxies, clusters, 
etc.) in the comoving volume dV at a comoving distance r from a given object A (quasar, 




(1) 
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galaxy, cluster, etc.). We concentrate ourselves on the plane perpendicular to the line of 
sight at the position of object A with redshift z d . The expected number of objects B within 
a radius d( of C from A is 

[l + Z AB (r)(l + z d y}-^==z, (2) 

where ubo is the present number density of objects B and e accounts for the evolution of the 
correlation function. In the following we consider only those excess population of objects 
B relative to the "background" ones, i.e., we exclude the contribution of the mean number 
density of objects B. 

If we assume that the mass density profile of object B has the form of a singular 
isothermal sphere (SIS) with velocity dispersion erg, which is a good approximation for the 
dark matter distribution in galaxies and in clusters of galaxies, the mass contribution from 
B at position ( to an area of tiQ is simply 

^Co ft~ 1[Co/0 <C,0)d9, Co<C; 



m B (CoX) = < 



(3) 

Co £ <C,0)d9, Co > C, 



in which k((, 9) = \Jl — (C/Co) 2 sin 2 9. The expected mass contribution to the area nQ 
from an ensemble of objects B is thus 

m(Co) = / m B (Co,C)^(C)- (4) 

Therefore, we can obtain the expected mean surface mass density of objects B around a 
given object A simply by m(£o)/7rCo% which reads 

WCo) = 4^o(l + ztfrl (^j F(Co), (5) 

where 

rCo fir rr max /-sin -1 (Co /C) 



/■c.0 fir rrmax /-sin ^(,o/tJ 

/ x(CMC/Co) / «(C, o)M + 2 / x(CKC/Co) / <C,e)d9 

Jo Jo Ji Jo 



(6) 
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and 

/ j. i (C/Co) 7 A T 7-i,/ r 2_i ' r min < C — r maxi 

x(Co, C) = < rmoB/( ^ (7) 

[ (C/Co) 7 JrZinK C < r mm- 

Suppose that the objects 5 follow a luminosity distribution function <Pb(Lb)cILb and 
the luminosity L B is related to the velocity dispersion through 

in which L* B and are the characteristic luminosity and corresponding velocity dispersion, 
respectively. The total surface mass density given by £ab(?") is the sum of all the objects B 
over their distribution of velocity dispersion or luminosity: 

Sab(Co)=4(1 + ^)V 2 ^F(Co) j / o (^J MLs)dL B . (9) 

However, applying SIS model to an arbitrarily large radius may lead to an overestimate 
of the matter contribution of objects B. A reasonable hypothesis is that the SIS profile is 
truncated at a radius Rb so that the total mass of object B is M B = 2R B cr B /G- So, the 
mass contribution m B (CoX) [ eo ±-(3)] of object B at C to an area of irQ should be properly 
replaced by uib^Rbi Co, 0- This gives a mean surface mass density of objects B around the 
object A to be 



Sab(Co) = 4(1 + ^)V 2 ^F(Co) jf° i^j ML B )dL B , (10) 



where 



F{(o) = 2— — / x Co, Qd C/Co)- (11) 

Co \ r 0/ J M B 

If we introduce the mean spatial mass density of objects B, Q B , in unit of the critical mass 
density po of the universe, eq.(10) can be written as 

Z AB ((o) = 4fi B po(l + *d)V /(Co), (12) 
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in which /(Co) = (r /2R B )F(( ). This formula was used by Wu & Fang (1996) for 
the estimate of the matter contribution to the cluster environments from the cluster- 
cluster correlation function. The advantage of this expression is that the effects of the 
observationally determined luminosity function (p B (L B )dL B and of the experiential formula 
eq.(8) are represented by a single parameter fl B . 

It appears that -F(Co) or /(Co) is a slowly varying function of Co if 7 is close to 2. As 
a result, the surface mass density around a given object provided by an ensemble of SIS 
objects that obey a correlation function of £(r) ~ r~ 2 is nearly a constant, which will be 
numerically verified below for various objects. The simple reasons are as follows: The total 
population of objects B on the plane perpendicular to the line of sight goes as N B ~ C if 
their spatial two-point correlation function has form of £(r) ~ r -2 . While the projected 
mass tub of a SIS is proportional to C, the expectation of the mass contribution by all the 
objects B is then m ~ m B NB ~ C 2 - Finally, the expected surface mass density of objects B 
around A, T^ab = m/ 7rC 2 , remains roughly unchanged. 

We now have two ways to evaluate the surface mass density around object A 
contributed by the two-point correlation function £ab(?~): (I)The conservative estimate 
of Tjab [eq.(10)], provided by a SIS density profile for B with a proper radius cutoff and 
an observationally determined luminosity function for the distribution of population B; 
(II)The optimal estimate of Y, A b [eq.(12)], assuming again a truncated SIS density profile 
for B but characterizing the distribution of all the population B by their matter density 
parameter fl B in the universe, especially when Q B = Q where fl is the present mean 
mass density parameter of the universe. We will employ these two models to compute the 
matter contributions from galaxy-galaxy, cluster-cluster, cluster-galaxy and quasar-galaxy 
correlations. The cutoff radii of galaxies {R g ) and clusters of galaxies (R c ) are assumed to 
be R g « 0.2 hr l Mpc and R c « 1.5 h~ l Mpc, respectively, according to the suggestions 
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by dynamical analysis and observations (e.g. Bahcall, Lubin & Dorman 1995). In the 
numerical calculations below we further require that the minimum separation r min in the 
correlation function £ab(?") be greater than the cutoff radius of object A. 



2.2. Galaxy-galaxy correlation 

We take the two-point galaxy-galaxy correlation function £ gg (r) from the CfA survey 
(Davis & Peebles, 1983), which provides r = 5.4 ± 0.3 hr l Mpc and 7 = 1.77 ± 0.04 at 
scale 10 kpc< r < 10 h' 1 Mpc. We classify the galaxy population as the early-type galaxies 
E/S0 and the spiral galaxies S in terms of their morphologies and adopt the composition 
given by Fukugita & Turner (1991): (71,72,73) = (12 ± 2%, 19 ± 4%, 69 ± 4%), where 
{i} = (1,2,3) =(E,S0,S). The distribution of different galaxies with luminosity is described 
by the Schechter luminosity function 

4> g ,i{L 9 )dL g = <t>* g {L g /Lir a eM-L 9 /L* g:i )dL g /L* g:i , (13) 
where 0* = (1.56 ± 0.34) x 10~ 2 h 3 Mpc~ 3 , a g = -1.07 ± 0.05 and 
(<J* 2 , cr* 3 ) = (225l2o> 206l2o> 144 ±i3) km s" 1 - Alternatively, the index in eq.(8) 
reads v = 4 for E/S0 galaxies and v = 2.6 for S galaxies, respectively. 

Numerical computation shows that /(Co) ~ 8.6 and varies only 9% over distance 
Co = 0.005 - 0.5 hr x Mpc. Therefore, the overall surface mass density S 99 (Co) from the 
galaxy-galaxy correlation around a given galaxy can be approximately taken to be a 
constant matter sheet: 

{ 0.0027^ (1 + h (F/0.64) g cm" 2 , (I); 
^99 = \ (14) 
[ 0.011Q fl (1 + z d y h (f/8.6) g cm" 2 , (II), 

in which the total error includes the errors in the adopted galaxy luminosity function and 
the galaxy-galaxy correlation function. Combining the above results obtained from the two 
models yields il g 0.25, in good agreement with the claim by Bahcall et al. (1995). 
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2.3. Cluster-cluster correlation 

Cluster-cluster correlation function £ cc (r) is represented by r = (20 ± 4.3) hr x Mpc 
and 7 = 1.8 over scale of ~ 5 hr x Mpc — ~ 75 hr x Mpc (Postman, Huchra & Geller 1992). 
We employ both the X-ray luminosity function and the mass function of galaxy clusters to 
estimate the cluster matter contribution from £ cc (r) to a given cluster at Zd- 

The X-ray luminosity function of clusters of galaxies shows a strong evolution with the 
cosmic epoch. Within the low redshift of z d ~ 0.15 (Edge et al. 1990) 

<f> c (L x )dL x = <j)* c L~*l exp(-L xM / L* xM )dL xM , (15) 

in which 0* = i()- 6 - 57±012 h 3 50 Mpc" 3 , a c = 1.65 ±0.26, L* ;44 = B.ltH and L xM is the X-ray 
luminosity in units of L 44 = 10 44 erg s -1 . For the distant clusters of galaxies at redshift 
ranging from 0.14 to 0.6 (Henry et al. 1992), 

(j) c (L x )dL x = <p* c L x * A ldL xM . (16) 

For the three redshift shells {(0.14, 0.20), (0.20, 0.30), (0.30, 0.60)}, a c = 
{2.19±0.21, 2.67±0.26, 3.27±0.29} and 0* = {5.85±0.25, 6.82±0.51, 12.33±3.87} x 10" 7 h 3 50 
Mpc _3 (L 44 ) Qc_1 . We utilize the quasi-Faber- Jackson relation for the X-ray selected galaxy 
clusters to convert X-ray luminosity into velocity dispersion (Quintana & Melnick 1982; see 
also Wu & Hammer 1993), 

L. = l<P- n ol" kg ergs- 1 . (17) 

The more recent optical and X-ray observations provide essentially a similar result (see Wu 
& Mao 1995 and references therein): L x = 10 32 64 a 4 erg s . Moreover, we introduce a 
cluster velocity dispersion cutoff a c>m i n at the faint end of the luminosity distribution to 
ensure the convergence of the integration / al4> c (L x )dL x . We choose <J c ,min = 508 km s _1 
which corresponds to an M* cluster (see below). The resulting surface mass density from 
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£ cc (r) in terms of model I reads 



S nn — < 



< z d < 0.15; 
0.14 <z d < 0.20; 
0.20 <z d < 0.30; 
0.30 <z d < 0.60, 



(18) 



0.66lg;ijjj x 10- 4 h (F/0.27) (1 + z d y g cm- 2 , 
2.0lJ:| x 10- 4 h (F/0.27) (1 + z d f g cm" 2 , 
2.81^ x 10~ 4 h (F/0.27) (1 + z d f g cm" 2 , 
8.3lt;l x 10~ 4 h (F/0.27) (1 + z d y g cm" 2 , 

in which the variation of / = (r /2R c )F is less than 10% from / = 1.8 for ( = 0.5 - 10 hr 1 
Mpc. 

A relatively simple way of estimating E cc is to use the cluster mass function established 
by Bahcall & Cen (1993): 



n c (> M c ) = <P* C (Mc/M*)" 1 exp(-M c /M c *), 



(19) 



where 0* = 4 x 10~ 5 h 3 Mpc -3 and M* = 1.8 x 10 14 hr 1 M Q . M c refers to the cluster mass 
within R c — 1.5 hr x Mpc radius sphere of the cluster center, which is related to our cluster 
model through M c = 2a^R c /G. A straightforward computation gives 



(20) 



in which we have utilized the limit a c > a c min 



1-318:1 x 10- 4 (1 + z d Y h (F/0.27) g cm" 2 , (I); 
0.0083fi c (1 + z d y h (f/1.8) g cm" 2 , (II), 

We have also shown the result of model II, 
represented by the cluster matter density parameter fl c in the universe, which is consistent 
with the finding of Wu & Fang (1996). Apparently, the result given by model I in the above 
equation is compatible with eq.(18) derived from the cluster X-ray luminosity function at 
low redshift z d < 0.3. Note that a comparison of the results by the two models in eq.(20) 
yields a relatively small value of the cluster matter density: VL C = 0.016. This is primarily 
due to the adopted low limit a c > cr c min = 508 km s _1 . We will not explore how <7 c , m in 
affects the estimate of E cc because the maximum cluster matter contribution can be figured 
out simply by setting Q c = Q in eq.(20). 



- 12 - 



2.4. Cluster- galaxy correlation 

It has been found that the shape of the cross correlation function of Abell clusters with 
Lick galaxies (, cg ( r ) is slightly steeper than i gg {j-) and £, cc ( r ), with 7 ranging from 1.7 to 2.5 
(Lilje & Efstathiou, 1988), while the amplitude of £ C9 (r) has not been well constrained. We 
adopt the approximate form by Peebles (1993): r = 15 ± 3 hr l Mpc and 7 = 2, which is 
available on scale roughly 0.5 hr l Mpc< r < 40 hr l Mpc. We compute both the matter 
contribution S ffC of clusters surrounding a given galaxy and the matter contribution S C9 of 
galaxies surrounding a given cluster. Replacing £ gg in eq.(14) and £ cc in eq.(20) by £ C9 gives 

S sc = 

and 

Here / appears to be nearly a constant over scale 0.01 h~ l Mpc< ( < 10 hr l Mpc 
in eq.(21). However, because of the small radius cutoff of density profile and the r~ 2 
distribution of spatial number density, galaxies lead to a relatively large variation of up to 
48% from / = 3.7 for ( ranging from 0.5 to 10 h~ l Mpc in eq.(22). 

2.5. Quasar-galaxy correlation 

Observations have shown that quasars are preferentially located in region of higher 
than average galaxies, usually in small groups of galaxies, and behave in a similar way to 
the galaxy-galaxy correlation: i qg {r) ~ £ gg (r) (Bahcall & Chokshi 1991; references therein). 
This arises from the fact that both quasars and galaxies trace the same large-scale structure 
of the universe. The proportional coefficient up to r = 0.25 h^ 1 Mpc is roughly 2.3 for the 
optically-selected quasars at z d < 0.5, and 2.8 and 8 for radio-loud quasars at z d < 0.5 



2.7tl\l x ID" 4 (1 + z d y h (F/l) g cm" 2 , (I); 
0.017Q c (l + z d yh(f/5)gcm- 2 , (II), 



0.0031^ (1 + ^)^(^/0.1) gem- 2 , (I); 
0.013^ (1 + z d y h (f/3.7) g cm" 2 , (II). 



(22) 
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and z d ~ 0.6, respectively. Nevertheless, we assume that £ 9ff (r) continues as ~ r~ L8 to a 
maximum separation of 10 h~ l Mpc as for the galaxy-galaxy correlation. We have checked 
an alternative model proposed by Bahcall & Chokshi (1991) that (, qg (r) has a steeper 
shape with 7 = 2.5 beyond and also normalized at r = 0.25 hr l Mpc. It is found that 
these two assumptions have resulted in the roughly same value of /. Nevertheless, the 
minimum separation between the quasar and its neighbor galaxies can be in principle taken 
to be r min ~ 0. This leads to / = 39 with a variation of smaller than 28% over distance 
0.01 < Co < 0.1 h^ 1 Mpc. Thus, the galactic matter contribution around a quasar can be 
obtained simply by increasing T, gg by a constant proportional factor: 



2->qg ~ 



2.3rE g9 , for optical quasars at z d < 0.5; 

2.8rE 9g , for radio loud quasars at z d < 0.5; (23) 
8rE S9 , for radio loud quasars at Zd ~ 0.6, 



where E S9 is given in eq.(14) and r accounts for the change of the function F or /: 
r = 39/8.6. The very recent observations of 20 luminous quasars at z < 0.3 with HST 
(Fisher et al. 1996) have found that the ratio of t; qg to £ gg is 3.8 ± 0.8. This factor is 
basically consistent with the result of Bahcall & Chokshi (1991). 



3. Applications in Statistical Lensing 

3.1. Gravitational lensing by a uniform matter sheet 

A uniform matter sheet E was introduced by Turner, Ostriker & Gott (1984) to model 
the cluster matter surrounding a galaxy. Basically, gravitational lensing by a massive object 
having a deflection angle oc with assistance of a uniform mass sheet E is described by the 
following lensing equation: 

f 1 " y^V " = TT«. ( 24 ) 
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where 

c 2 D s 

KrU = ^GD^D, (25) 
is the critical surface mass density, (3 and 8 are the positions of the source and of the 
corresponding images, and D d , D s and D ds are the angular diameter distances to the 
lens, to the source and from the lens to the source, respectively. It appears that the 
contribution of a matter sheet to gravitational lensing becomes to be significant only when 
£ reaches a value comparable to S crii . S crii tends to infinity at both D d = [z d = 0) 
and D d = D s (z d = z s ). For a typical source (e.g quasar) at z s = 2 the minimum value 
of S cr j t is 0.82 h g cm~ 2 . It turns out, according to the calculations made in the above 
section, that the clustering of galaxies and galaxy clusters can provide a matter density 
of as high as S/E crit ~ 1CT 2 around a given galaxy or cluster, while the galactic matter 
around a given quasar can reach E gs /S cr jt ~ 1CT 1 , provided that Vl g ~ f2 c ~ f2o = 1- 
This environmental matter may raise the amplitude of gravitational lensing associated 
with foreground galaxies, clusters of galaxies and quasars. In the following subsections we 
explore one of the consequences arising from magnification bias due to the presence of the 
clustering of galaxies and clusters of galaxies: the associations of the angular positions of 
background quasars with foreground galaxies, clusters of galaxies and quasars. To do this 
we will adopt a stable clustering model e = —1.2, which has been shown to be consistent 
with the universal relation ^As( r ) ~ (r/r ) -1 ' 8 for clustering of galaxies, clusters of galaxies 
and quasars (Bahcall & Chokshi 1991). Alternatively, the surface mass density given by 
£,AB{r) refers to the value over a comoving area. A factor of (1 + z d ) 2 should be multiplied 
to get the physical surface mass density in an Q = 1 universe. 



- 15 - 



3.2. Quasar enhancement factor 

Quasar overdensity behind a massive object is usually characterized by the quasar 
enhancement factor q, which is the ratio of the observed quasar surface number density over 
the association area around the foreground object to the intrinsic one N(< m b ) above a 
limiting magnitude m b . In practice, N(< m b ) is obtained using the quasar number density 
in the remaining area far away from the foreground object. Moreover, the "unaffected 
background hypothesis" is often employed, i.e., the observed quasar number count is not 
significantly different from the intrinsic one N(< mi,). In order to compare with the 
observations that search for quasars over an area with distance ranging from 6\ to 62 from 
the foreground object, we use the average enhancement factor q: 

q = w^wr (26) 



where (Narayan 1989) 

N(< m b + 2.5 log fj,) 1 



(27) 



N(<m b ) /i' 

and \x is the lensing magnification introduced by all the gravitational matter associated with 
foreground object. Eq.(27) has accounted for both the magnification bias (2.5 log /x) and 
the area distortion (l//i) because of light deflection near the lens. For the optically-selected 
quasars we adopt the number-magnitude relation from Boyle, Shanks & Peterson (1988) 

N(< m h ) = 4.66 x io°- 86 ( m ^ 19 - 15 ) < 19 15 

K J (28) 
N(< m b ) = -10.95 + 15.61 x io - 28 ^- 19 - 15 ), m b > 19.15, 

which is valid within m b < 21 and z < 2.2. If the foreground object has a spherical matter 
distribution, the magnification induced by this lens with an additional uniform mass sheet 
E becomes 



r — \ 

1 - 



S V I ^ _ 0,: <>(0) \ L _ 0e da(6) \ 



^critj \ 6 ol(6 e )) V ot(9 E ) d6 J 
where Qe is the Einstein radius corresponding to (3 = in eq.(24). 



(29) 
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3.3. Quasar-galaxy associations 

Statistical evidences for quasar-galaxy associations have been accumulated for a decade 
since the first observation by Tyson (1986). It is generally believed that statistical lensing 
by galaxies as well as their host clusters is most likely to be the cause. Here, we investigate 
how galaxy environments, namely, the galaxy-galaxy correlation and the cluster-galaxy 
correlation, affect the lensing properties of quasar-galaxy associations. We again utilize a 
SIS model for the mass density profile of a galaxy and the same morphological composition 
7i (i = 1,2,3) for galaxy population as we adopted in section 2.2. For the galactic matter 
(SIS) with a uniform matter sheet S 9 = T, gg + T, gc provided by i gg (r) and £ C9 (r), eq.(29) 
reduces to 

"=<^rrk^ (30) 

where 6e = ^b(1 — S 9 /S cr j t ) _1 and 6 9 e = 4n((j g /c) 2 (D ds / D s ) are the the Einstein radius 
with and without E g , respectively. The expected quasar enhancement factor (q) is obtained 
by averaging q over the luminosity and spatial distributions of galaxies: 

So' ^D 2 d (l + z d f(dr proPtZ Jdz d )dz d l^f q { 'y i <j> i {L g )dL g J 

(?) = 7^= f . (3D 

Io s ^D 2 d (l + z d f(dr proP)Z J dz d )dz d ^ / ^{L^dLg 

in which dr proPjZd is the proper distance within dz d of z d , L min is the luminosity of the 
faintest galaxies in the sample and a constant comoving number density of galaxies has 
been presumed. In effect, the denominator is the total number of galaxies with L > L min 
and z < z s . Note that the observations actually measure the surface number density 
of low-redshift galaxies around high-redshift quasars, giving rise to the so-called galaxy 
enhancement factor (Wu 1994). In principle, the foreground galaxy enhancement factor 
should be equal to the background quasar enhancement factor q in the sense that one 
counts only the quasar-galaxy pairs within a given distance range in a complete sample. 
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A detailed numerical computation of eq.(31) using the parameters in the observations of 
quasar-galaxy associations has been recently made by Wu, et al. (1996). We present in 
Table 1 the observed and expected quasar enhancement factors behind galaxies with and 
without the environmental matter H g . It appears that, as compared with galaxies alone as 
lenses, the mean quasar enhancement factor (q) is typically increased by ~ 0.01 and ~ 0.1 
for model I and model II (assuming that Q g = Q c = Qq = 1), respectively. Therefore, (q) 
is little affected by the galaxy environments if galaxies and clusters of galaxies contribute 
no more than < 25% to the matter of the universe (model I), while the theoretical 
predictions fit better the measurements of the quasar-galaxy associations with the help of 
the environmental matter if Q g ~ Q c ~ Q = 1 (model II). 

EDITOR: PLACE TABLE | HERE. 

3.4. Quasar-cluster associations 

A significant quasar overdensity has been recently detected on scale of ~ 10 arcminutes 
behind foreground clusters (Rodrigues- Williams & Hogan 1994; Wu & Han 1995; 
Rodrigues- Williams & Hawkins 1995; Seitz & Schneider 1995). It has become clear today 
that cluster matter alone cannot account for the large quasar enhancement in the scenario 
of gravitational lensing, and large-scale matter clumps traced by clusters of galaxies should 
be invoked. Wu & Fang (1996) have analyzed the four measurements of the quasar-cluster 
associations and found that the surface mass density required to produce the reported 
quasar enhancement factors should be as high as ~ 0.2 h g cm -2 . They only computed the 
matter contribution S cc from the cluster-cluster correlation function £ cc (r) and concluded 
that the resulting E cc is at least an order of magnitude lower than ~ 0.2 h g cm -2 even if 
the total cluster matter is large enough to close the universe. This has been clearly shown 
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in eqs.(18) and (20). The present paper has further demonstrated that the discrepancy 
cannot be resolved by involving the cluster-galaxy correlation £ C9 (r): the galaxies in the 
proximity of clusters of galaxies following £ C9 (r) provide a mean surface mass density S C9 
[eq.(22)] comparable with E cc [eq.(20)], which is insufficient to explain the associations of 
the distant quasars at z s ~ 1 - 2 with the nearby clusters at z d ~ 0.15. 



3.5. Quasar-quasar associations 

The presence of close quasar pairs with separation of a few arcseconds but with different 
redshifts has been a puzzle for over two decades. Four events have been so far detected 
among ~1000 - 2000 surveyed quasars (Burbidge et al. 1996). While the probability 
of finding four quasar pairs within 5" based on accidental configuration turns out to be 
relatively small, it is argued that quasars probably have significant non-cosmological redshift 
components. An alternative to this non standard model is to attribute the quasar-quasar 
associations to magnification bias by foreground matter inhomogeneities, as was suggested 
by Schneider (Burbidge et al. 1996). However, there have been no any realistic lensing 
models which could be constructed to produce the reported quasar-quasar associations. 
Now, the quasar-galaxy correlation function (, qg (r) may provide a solution. 

Quasars are preferentially located in the over-density regions which have a mean surface 
mass density T> qg . Thus, the probability of detecting a background companion quasar near 
a foreground one within a given range is a factor of q higher than that over the similar area 
in the rest of the sky, where the enhancement factor q is given by eq.(27). Note that q is 
independent of the magnitude of foreground quasar, which goes down to 20.7 in the four 
examples of the quasar-quasar associations. All the four background companion quasars 
in the reported events are brighter than ~ 19.15 in B magnitude, the turn-over point in 
the quasar number- magnitude relation [eq.(28)]. This removes the concern by Schneider 
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(Burbidge et al. 1996) that quasar counts are too flat at faint magnitude to provide a large 
enhancement factor. 

It is rather difficult today to figure out the statistical significance for the reported 
quasar-quasar associations since they did not result from a uniform quasar sample with 
a definite limiting magnitude. Furthermore, the quasar-galaxy correlation increases 
considerably with cosmic epoch, which has not been well constrained. While the foreground 
quasars in the four known examples have redshifts ranging from 0.44 to 1.62, we take the 
relation (, qg (r) ps 8£ S9 (r) established so far at the highest redshift Zd ~ 0.6 to estimate 
its contribution to gravitational magnification. Recall that the recent HST observations 
have not found the significant difference in the clustering amplitude between the radio-loud 
and radio-quiet quasars (Fisher et al. 1996). Our numerical computation shows that the 
expected q for the four quasar-quasar association events ranges from 1.2 to 1.5 in model 
I and from 2 to 8 in model II, where we have utilized the mean value of £ gg (r) and given 
no information about the uncertainties due to the ignorance of the errors in the adopted 
£, qg {r). It turns out that about 2-8 times more distant quasars would be found near the 
positions of foreground quasars than in the rest of the sky, provided that Q g ~ Q = 1 
(model II). This enhancement factor essentially accounts for the probability of detecting 
four quasar companions within 5" of the primary quasars if the surveyed sample consists 
of as large as 2000 quasars. Recall that the expected quasar pairs based on the random 
distribution is about 0.6 (Burbidge et al. 1996). However, if £l g ps 0.25, as indicated by 
current observations, it is unlikely that the situation of the quasar-quasar associations would 
be significantly improved by introducing the lensing effect due to the quasar environmental 
matter. 
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4. Discussion and conclusions 

The connection of galaxies, clusters of galaxies and quasars with their environmental 
matter is characterized by their spatial two-point correlation function £(r). Assuming a 
SIS model for matter distribution in galaxy and cluster of galaxies, we have computed the 
expected matter contributions of the surrounding galaxies and galaxy clusters following £ (r) 
to a given galaxy, cluster or quasar. We have then employed these surface mass densities 
from £(r) for the estimate of the background quasar enhancement factor in the vicinities of 
foreground galaxies, clusters and quasars according to gravitational lensing magnification 
bias. 

Significance of the environmental effect in study of gravitational lensing is demonstrated 
by the ratio of the surface mass density £ provided by £(r) to a critical mass density £ crit . 
It appears that the environmental galaxy/cluster matter surrounding a galaxy /cluster 
has S/S cr j t < 10~ 3 , showing a negligible contribution to gravitational lensing, if their 
mass densities are no more than ~ 25% of the mass density of the universe (model I). 
In this circumstance, it is unlikely that the reported quasar-galaxy, quasar-cluster and 
quasar-quasar associations can be attributed or fully attributed to gravitational lensing 
by the environmental matter. However, in the extreme case of Q g ~ Qq — 1 and/or 
fi c ~ fi = 1, the amplitudes of the quasar-galaxy and quasar-quasar associations can be 
increased by their neighbor galaxies and clusters of galaxies, resulting in an expected value 
that is basically in agreement with the observed quasar overdensity behind foreground 
galaxies and quasars. Nevertheless, in no case can the lensing magnification by clusters and 
their surrounding matter account for the recently detected quasar-cluster associations on 
scale of ~ 10 arcminutes, in consistent with the claim of Wu & Fang (1996). 

It is noted that the two-point correlation function £(r) is observationally and 
statistically established. Therefore, our conclusions hold true only for an ensemble of 
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objects. It is just appropriate to the associations of background quasars with foreground 
galaxies and clusters of galaxies. However, the quasar-quasar associations might not 
result from a complete sample and thereby, our estimate about the background quasar 
enhancement factor near foreground quasars should be employed cautiously for the present 
data. 

It should also be noticed that the critical surface mass density S cr jt depends on the 
cosmological models. We have only adopted a matter dominated flat universe in the 
present paper. If the cosmological constant Ao is nonzero, as was suggested recently by 
the measurement of the Hubble constant H ~ 80 km s _1 Mpc -1 (Ostriker & Steinhardt 
1995; references therein), gravitational lensing by £(r) would be more efficient than it is in 
a A = universe. For instance, the minimum S cr j t reduces to 0.57 h g cm -2 for z s = 2 
in an Q Q = 0.2 and A = 0.8 universe. It deserves to be investigated how the cosmological 
constant affects the present calculations. Moreover, S crii appears to be smaller in a low 
density universe. So, if Qo ~ Q g = 0.25, one would expect a larger environmental effect 
than what we have discussed here for a flat cosmological model. 

Our present results may suffer from uncertainties in the evolutionary models of the 
correlation function with cosmic epoch. Observations have indicated a significant evolution 
of galaxy-galaxy correlation: Both the correlation length r and the index e of eq.(2) are 
found to change with redshift (Shepherd et al. 1996; reference therein). Recall that for 
quasars at z s ~ 2 the lensing galaxies are most likely to be at z d ~ 0.5, where the evolution 
effect cannot be neglected. Alternatively, the evolution of the quasar-galaxy correlation 
has been seen in eq.(23), while foreground quasars in the study of the quasar-quasar 
associations are actually beyond these redshift ranges. Therefore, our predicted amplitude 
of the quasar-quasar associations may have large errors. Fortunately, the quasar-cluster 
associations are little affected by the evolution since all the clusters in the measurements of 
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the quasar-cluster associations are at relatively low redshift z d ~ 0.15. It may be interesting 
to study the dependence of the quasar-galaxy, especially quasar-quasar associations on the 
matter clustering at moderate/high redshift, which may offer a possible way to test the 
evolution models of the formation of large-scale structions. 

Finally, a thin-lensing hypothesis has been used throughout this paper, which is 
apparently not a good approximation for large-matter clumps in the universe. We should 
thus point out that the results obtained from a simple analysis of the two-point correlation 
functions provide only a useful estimate of the amplitude of gravitational lensing generated 
by large-scale structures on various scales. An extensive study with N-body simulations 
may give a new sight into the subject. 

We thank an anonymous referee for valuable comments and suggestions. WXP was 
supported by the National Science Foundation of China and a World Laboratory fellowship. 
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Table 1: Quasar-galaxy associations. 



observations* expectations 



authors 


lobs 


Qgalaxies 


9(1) 


q(ll) 


Crampton 


1.4 ±0.5 


i 1 i +0.02 
i - ii -0.04 


i i O+0.03 
i - iz -0.04 


i 9 o+0.03 


Kedziora-Chudczer 


~ 1 


i ni +0 - 00 

i - ui -0.00 


i.0218:S 


i i i +0.01 
i - ii -0.00 


Magain 


~ 2.8 




i.49l8:l^ 


1-6418:11 


Thomas 


1.7 ±0.4 


i-08±8:g 


1-0918:81 


1.20+8:81 


Van Drom 


~ 1.46 


1 04+ - 01 

-■-•^-0.03 


i - ud -0.02 


-L-J-O-O.Ol 


Webster 


~ 2 


i - uo -0.02 


-L-U' -0.02 


i i O+0.02 
1 - 1 ^-0.02 




1.0 ±0.3 


l.UO_ 02 


i - uo -0.02 


i.iol8:81 


Yee 


1.0 ±0.2 


i n3 +0 - 01 


1 04+ - 01 

-L-U^-O.Ol 


1.0918:81 




0.9 ±0.1 


1 02+ - 01 
-■-•u^-o.oo 


1 ns +0 - 01 

i - u,3 -0.01 


1-0818:85 



*Data are taken from Narayan (1992) and Wu et al. (1996). 
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